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P-modes in rapidly rotating stars - 

looking for regular patterns in synthetic asymptotic spectra 
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According to a recent ray-based asymptotic theory, the high-frequency p-mode spectrum of rapidly rotating stars is a 
superposition of frequency subsets associated with dynamically independent regions of the ray-dynamics phase space. 
At high rotation rates corresponding to typical 8 Scuti stars, two frequency subsets are expected to be visible : a regular 
frequency subset described by a Tassoul like formula and an irregular frequency subset with specific statistical properties. 
In this paper, we investigate whether the regular patterns can be detected in the resulting spectrum. We compute the 
autocorrelation function of synthetic spectra where the frequencies follow the asymptotic theory, the relative amplitudes 
are simply given by the modes' disk-averaging factors, and the frequency resolution is that of a CoRoT long run. Our first 
results are that (i) the detection of regular patterns strongly depends on the ratio of regular over irregular modes, (ii) low 
inclination angle configurations are more favorable than near equator-on configurations (iii) in the absence of differential 
rotation, the 2Q. rotational splitting between m — 1 and m = — 1 modes is an easy feature to detect. 
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1 Introduction 

The launching of the space missions COROT and KEPLER 
is bringing a new wealth of observational data on the as- 
teroseismology of many different types of stars. For slowly 
rotating stars, e.g. the sun, the approximate spherical sym- 
metry of the system enables us to classify modes accord- 
ing to well-defined sets of integers. This mode structure al- 
lowed us to extract information on the stellar interior from 



dependent subspectra, corresponding to different phase space| 
zones for the ray dynamics. The subspectra corresponding 
to integrable zones give rise to regular sequences of frequen- 
cies, whereas a superimposed subspectrum is associated to 
chaotic dynamics and displays regularity only in a statistical 



the sp ectrum of oscillation modes dChristensen-Dalsgaard , 
2002) . 



However, for rapidly rotating stars, the deformation of 
the star breaks the spherical symmetry. It has recently been 
shown that perturbation theory fails to yield accurate pre- 
dictions of osc illation spectra even for stars with mod erate 
rotation rates dLignieres et all 120061 Reese et al. , 20061) . In 
Lignieres & Georgeot! d2008[ 120091) . an asymptotic theory 
of spectra of such stars was built, based on acoustic ray dy- 
namics. It was shown that this dynamics is Hamiltonian, and 
undergoes a gradual transition from integrability to chaos 
when the rotation increases. At moderately rapid rotation, 
the phase space is divided into integrable and chaotic zones, 
w here the dynamics is qualitatively d ifferent. It was shown 



In order to connect such results to observed spectra, sev- 
eral questions have to be explored. The first one is the va- 
lidity of the asymptotic theory to the finite range of spectra 
w hich can be observed. A first answer was given in 



Lignie res & Georgeo t (2008, 2009), where it was shown thafl 



in 



Lignieres & George ot (2008, 2009) that this structure mod 



ifies the structure of the oscillation spectra. The acoustic 
mode frequencies cannot be in general associated to well- 
defined sets of integers. Instead, the spectrum is split into in- 



* Corresponding author: e-mail: francois.lignieres@ast.obs-mip.fr 



the predicted structure can be found in numerically 
computed relatively low-frequency p-modes of polytropic 
stellar models. An other important question is to separate 
these different subspectra in observational data, where the 
visibility of the modes plays an important role. A first step 
in this direction is to construct synthetic spectra based on 
-the asymptotic theory where visibility varies from mode to 
mode and try to extract from them important features of the 
system, as if they were observed spectra. This is the strategy 
we follow here. We also note that a first evidence of regular 
patterns in a p-mode spectrum of a rapidly rotating star has 
been recently obtai ned by analyzing the spec tra of a 6 Scuti 
CoRoT target dGarcia He rnandez et al. . 20091) . 
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2 Construction of synthetic spectra in the 
asymptotic regime 

We have built synthetic spectra with the following formulas: 

S(v,i) = y i L(v,v j ,A j {i),T) (1) 



3 



L(u, v^Ajii), T) = Aj(i) exp [-(u - utf/T 3 



(2) 



where Uj are the frequencies of the modes and Aj(i) their 
amplitudes, which are functions of the inclination angle i. 
We have represented modes with Gaussian functions of 
width r to take into account the finite resolution of the spec- 
trum. We have considered in this work a resolution V corre- 
sponding to 150-day-long runs. Our main assumptions are 
that the f requencies follow the asymp totic behaviour de- 
scribed in lLignieres & George ot (2009) and that the relative 
amplitudes only depend on the relative mode visibilitites, 
the visibility being approximated by the disk-averaging fac- 
tor computed for high-frequenc y p-modes in a polytropic 



3 Autocorrelation of the synthetic spectra 

The autocorrelation of the synthetic spectra has been inves- 
tigated, for £1 = Q.QHk, by varying three parameters: the 
frequency range of the spectrum, the inclination angle of 
the star, and an amplitude threshold that only retains the 
the highest amplitude peaks of the spectrum. The parameter 
domain that we explored is the following : The frequency 
range spans ng large separations S n = 2A^, where ng has 
been varied from 1 to 6. The amplitude threshold is chosen 
to keep Ns frequency peaks per large separation interval; 
Ng has been varied between 10 and 100. Once ns and Ng 
are fixed, the total number of frequencies of the spectrum is 
ns x Ns- The inclination angle has been varied from to 
90°. 



stellar model (see formula (37) in 

Above a certain rotation rate, the p-modes spectrum is dom- 
inated by two families of modes, the island 
modes and the chaotic modes, as the visibilities of the other 
modes become negligible. The island modes are associated 
with island chain structures of the ray dynamics p hase space 
and th eir spectrum follows a simple formula dReese et all 



2008). 



lira 



= A„h + Ag£ + A m |m| — mCl + cj re f (3) 

where I < 3 and Aft, A/, A m depend on the star's structure. 
The quantum numbers h and t are defined from the spa- 
tial distribution of the modes in a meridional plane. They 
correspond respectively to the number of nodes along and 
across the stable periodic trajectory associated with the is- 
land chain. By contrast, chaotic modes are associated with 
chaotic regions of phase space, and their frequency spec- 
trum is said to be irregular because it is not described by 
a smooth function of 3 integers. The frequency spectrum 
has nevertheless specific statistical properties. Indeed, if we 
consider chaotic modes of the same symmetry class, the dis- 
tribution of the consecutive frequency spacings <jj = u>i+i — 
uji (scaled by the mean frequency spacing < cJi+i — u>i >) 
is close to the parameter- free Wigner distribution P(a) = 
7rcr/2 exp(— 7rcr 2 /4) in accordance with the prediction of 
Random Matrix Theory. A symmetry class corresponds to a 
given azimuthal number m and a given symmetry (+ or — ) 
with respect to the equator. For each symmetry class, a spec- 
trum is thus determined from a realization of the Wigner 
distribution and the chaotic spectrum is the superposition of 
all these spectra. Finally in order to construct the total spec- 
trum, we need to know the ratio between chaotic and island 
modes. This ratio tends to increase with rotation and, for the 
rotation rate considered in the follo wing 17 = 0.617 where 
Q,k = \/GM/R 3 , it is close to 3.7 dlignieres & Georgeot , 



2009). 



We find that the inclination angle plays an important role 
in the search for regular patterns because it affects the rel- 
ative visibility of the chaotic and island modes. Indeed, the 
disk-averaging factor of the chaotic modes tends to increase 
towards equator-on configurations while the disk-averaging 
Lignieres & Georgeotl(l200 9)lgtor of island modes tends to decrease. 

Let us first consider a high inclination configuration i = 
63°. In figures 1 and 2, the autocorrelation of the whole 
spectrum (top panel), the chaotic spectrum (middle) and the 
island spectrum (bottom) are displayed for Ns — 57 and for 
two different frequency ranges ns = 3 (Fig. 1) and ng = 6 
(Fig. 2). The autocorrelation is defined as in statistics: it is 
computed after substracting the mean and it is normalized 
by the variance. The autocorrelation is therefore comprised 
between - 1 and 1 . The vertical scale is different for the auto- 
correlation of the island spectrum because, as expected, the 
Tassoul like formula leads to strong peaks associated with 
the An regular spacing and smaller peaks corresponding to 
linear combinaisons of A^, A/, A m . 

No structure is seen in the chaotic mode spectrum ex- 
cept for a peaks at 217 (corresponding to the 20,/ 5 n = 0.836 
peak of Fig. 1). This lack of structure can be understood 
from the fact that the chaotic spectrum is a superposition 
of the subspectra corresponding to the different symmetry 
classes m . These subspectra are statistically independent 
and, as a consequence, the consecutive difference statistics 
of the chaotic spectrum is generally close to a Poisson statis- 
tics thus explaining the featureless autocorrelation. The 
pole-on configuration is an exception to this rule as only 
axisymmetric modes m = are visible in this case and the 
superposition of only two independent subspectra does not 
lead to Poisson statistics. The 217 peaks is due to a rota- 
tional splitting between m = 1 and m = — 1 modes. In the 
p-modes asymptotic regime, the Coriolis force has indeed 
a negligible effect on the frequency because its characteris- 
tic time-scal e is much larger the n the oscillation time-scale. 
As shown in lReese et al.l (120061) . this asymptotic property is 
actually already correct at relatively low frequencies. Since 
the centrifugal force does not distinguish between — m and 
m modes, their frequencies can be considered as degener- 
ate in the rotating frame thus leading to a 2to$7 splitting in 
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Fig. 1 Autocorrelation of a high-frequency p-modes 
spectrum resulting from the superposition of a chaotic mode 
spectrum and an island mode spectrum. The top panel dis- 
plays the autocorrelation function of the whole spectrum, 
while the middle and bottom panels shows the autocorre- 
lation of the chaotic and island mode spectra, respectively. 
The x-axis represents the frequency lag in units of the large 
separation S n . The inclination angle is equal to i — 63°, the 
frequency range corresponds to 3 large separations and the 
171 highest amplitude frequencies have been retained in the 
spectrum. The only significative feature of the autocorrela- 
tion function is the 251 peak also seen in the autocorrelation 
of the chaotic spectrum. 



the inertial frame. In addition, the number of m = 1 and 
m = — 1 modes in the spectrum is relatively high because 
their disk-averaging annulation effect is not strong (and less 
important than for higher m modes), thus leading to a sig- 
nificative peak at 2ft in the autocorrelation function. 

As shown in the top panel of Fig. 1, signatures of the 
island mode regular patterns are not apparent in the auto- 
correlation of the total spectrum. This is due to the fact that 
the ratio of chaotic modes over island modes is too high in 
this case. Increasing the frequency range to 6 large separa- 
tions (Fig. 2) enables us to detect a peak at the large sepa- 
ration 5 n (the peak seen at a scaled frequency lag of unity 
on Fig. 1 or Fig. 2). This is expected as a larger frequency 
range leads to the build up of the peaks associated with the 
periodic features of the frequency spectrum. For the present 
high inclination angle configuration, limiting the analysis to 
the highest amplitude peaks does not help because in aver- 
age the visibility of the island modes is not higher than the 
visibility of the chaotic modes. 

As illustrated in figures 3 and 4 for i = 30°, low inclina- 
tion configurations are more favorable in order to detect the 
island mode regular patterns. The autocorrelation functions 
are shown for a frequency range ng = 3 and for two dif- 
ferent amplitude threshold N s = 68 (Fig. 3) and N s = 18 
(Fig. 4). In both cases, features associated with the island 
mode regular patterns are detected in the total autocorre- 
lation function. This is because, as compared to the high 



Fig. 2 Same as figure 1 except that the frequency range 
of the spectrum now spans 6 large separations, the number 
of frequencies retained per large separation remaining the 
same. The autocorrelation spectrum now shows a peak cor- 
responding to the large separation. 
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Fig. 3 Same as Fig. 1. The inclination angle is now equal 
to i = 30°, the frequency range of the spectrum spans 3 
large separations and the 204 highest amplitude frequencies 
have been retained. The main correlation peak of the island 
spectrum can be retrieved in the autocorrelation of the total 
spectrum. 



inclination case, the proportion of island modes in the total 
spectrum has increased. The comparison of figures 3 and 4 
also shows that, contrary to the high inclination case, the is- 
land mode regular patterns are more easily detected if the 
analysis is limited to the highest amplitude modes. Finally 
we checked that, as expected, the signature of these patterns 
is stronger if the frequency range is increased. 



4 Discussion and conclusion 

Although based on strong simplifying assumptions, the 
present study should provide some clues to conduct a search 
for regular patterns in observed p-mode spectra of rapidly 
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gravity darkening effects. The mode inertia could also be 
taken into account in modelling the amplitudes. 
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Fig. 4 Same as Fig 3. except that only the 54 highest am- 
plitude frequencies have been retained. This helps to detect 
some of the correlation peaks that were not significant in the 
previous case. 
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rotating stars. For example, we find that limiting the search 
for regular patterns to the highest amplitude frequency peaks| 
is not necessarily helpful. This depends on the inclination 
angle. Generally, low inclination configurations are more fa- 
vorable because the island modes are relatively more visible 
than the chaotic modes. A robust feature of the autocorrela- 
tion functions computed in the present study is the 2£1 peak. 
This feature does not rely on the asymptotic regime assump- 
tion, it rather requires that the Coriolis force has a negligible 
effect which is already true at relatively low p-mode fre- 
quency. Possible effects of the differential rotation on this 
peak should nevertheles s be tested , for e xample using the 
formalism described in lReese et alJ (2009). 

The two most important assumptions in the construc- 
tion of the synthetic spectra concern the asymptotic regime 
and the mode amplitudes. We assumed that island mode fre- 
quencies strictly follow the asymptotic formula (01 while it 
is in fact only appro ximate due to finite waveleng th effects 
(see the dicussion in lLignieres & Georgeot! \2Q0% \ The de- 
partures from the asymptotic formula (01 have been deter- 
mined for different stellar models and fre quen cy ranges in 
Lignieres etail d2006l) . iReese et al.1 d2008h and lReese et al. 



(2009). As a first step, a simple way to model these non- 
asymptotic effects in the present study would be to decrease 
the frequency resolution of the synthetic spectra. A second 
step would be to co mpute numerically a complete spectrum 



using the code TOP dReese et all 12006). Such computations 



are nevertheless very time consuming. Thus, simple mod- 
els of the spectra like the one presented in this study will 
still be useful as they enable us to cover a wider range of 
parameters. In what concerns the amplitudes of the modes, 
we lack a consistent theory to describe their excitation and 
non-linear saturation in rapidly rotating stars. The present 
model could nevertheless be improved through better deter- 
minations of the mode visibilities, including more realistic 
stellar models as well as non-adiabatic, limb-darkening and 
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